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SENSITIVITY  STUDY  OF  RIFLE  GAS  SYSTEMS 


ABSTRACT 

Results  of  a  sensitivity  study  of  the  M-16  rifle  gas  system  are 
presented;  this  study  is  based  on  a  simulation  of  rifle  gas  system 
operation  developed  in  the  BRL.  The  calculations  indicate  that  thermo¬ 
dynamic  variables  in  the  bolt  carrier  cavity  are  only  weakly  sensitive 
to  variations  in  the  following  parameters:  i)  pressure  and  temperature 
in  the  gun  barrel  when  the  bullet  passes  the  port,  ii)  friction  in  the 
duct  flow,  and  iii)  frictional  resistance  to  motion  of  the  bolt  carrier. 
The  computational  results  are  sensitive,  however,  to  the  chosen  origin 
of  time  on  tne  oscillogram  showing  barrel  pressure  at  the  port  station. 

Graphs  are  presented  for  a  typical  round  illustrating  pressure, 
temperature,  density,  and  piston  motion  histories  for  M-16  and  AR-18 
rifle  gas  systems. 
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I .  INTRODUCTION 


It  is  desirable  to  simulate  the  operation  of  the  gas  system  which 
extracts  a  spent  cartridge  and  positions  the  next  round  in  fthe  chamber 
of  an  automatic  gas-operated  rifle.  The  operation  cycle  is  a  rather 
complicated  gasdynamic  process  and  must  necessarily  be  described  with„a 
considerably  simplified  model. 

The  model  developed  in  Reference  1*  has  proven  to  be  reasonably 
accurate  in  predicting  pressure  history  for  the  bolt  cavity  of  the  M-16 
rifle.  This  model  is  applicable  to  the  firing  of  all  rifles  during  the 
initial  interval  before  moving  parts  interact  in  the  reloading  mechanism. 
This  interval  covers  a  small  portion  (  -1/2 Oth)  of  the  total  operation 
time,  but  during  this  period  the  pressure  and  temperature  attain  their 
peaks  and  the  system  receives  all  its  momentum  from  the  powder  gas. 

The  above  model  forms  the  basis  for  this  report.  In  particular, 
we  shall  study  the  sensitivity  of  the  gas  system  to  variations  in  input 
data.  Since  it  is  not  possible  to  survey  in  detail  the  many  factors 
which  affect  the  performance,  attention  will  be  concentrated  on  several 
of  the  more  notable  effects  caused  either  by  change  of  conditions  inside 
of  the  gun  (e.g.,  friction  in  gas  duct,  change  of  amnunition)  or  un¬ 
certainties  resulting  from  measurements  (e.g.,  aero  time  and  peak  pressure 
at  gas  duct  entrance) . 

The  basic  analysis,  assvmtptions,  and  approximations  are  described 
thoroughly  in  Reference  1  and  will  not  be  repeated  in  detail  here.  A 
brief  outline  of  the  model  is  given  in  Section  II.  Results  of  com¬ 
putations  and  measurements  for  the  M-16  rifle,  including  data  for  the 
sensitivity  study,  are  presented  to  illustrate  general  behavior  of  gas 
systems.  Finally,  data  are  presented  for  the  AR-18  rifle,  and  some 
comparisons  with  M-16  results  are  made. 

•Current  confutations  eufloy  the  "area  discontinuity"  treatment  of  the 
gun  tube  port ,  treated  in  the  ADDENDUM  of  Reference  1.  References 
are  listed  on  page  3?. 
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II.  THEORETICAL  MODEL  OF  GAS  SYSTEM 

For  the  early  stage  of  the  operation  cycle,  rifle  gas  systems  are 

essentially  equivalent  to  the  system  represented  by  the  idealized 

geometry  shown  in  Figure  1.*  Hot,  high-pressure  powder  gases  are 

extracted  from  the  rifle  barrel  and  fed  through  a  duct  into  a  chamber, 

where  the  gases  expand,  thereby  accelerating  a  piston.  After  the  piston 

has  traveled  a  certain  distance,  vent  holes  are  laid  free  and  the  gas  in 

the  chamber  expands  to  atmosphere  .  The  pressure,  temperature  and  density 

of  the  powder  gas  in  tiw  barrel  —  p  (t) ,  0  (t) ,  p  (t)  —  are  known,  and 

9  9  9 

constitute  the  input  to  the  problem. 

The  quantities  to  be  determined  are  the  flow  variables  in  the 
chamber  and  the  motion  of  the  piston.  The  equations  to  be  solved  are:* 


A  p  u  »  (V  ,  +  Ajt  ) (dp/dt)  +  Apv  +  venting  term  (1) 

eeecxcB  ccb 

m  h 


-  .  .  -  d/dt  [pSY  +  1/2  S  v  2  V  +  A  P„v_dt]  + 

e  tot  ccc  c  B  c  c  o  c  B 


(2) 


venting  terms 


aVd‘  “  5c*c"b  "l  +  *  "b 


(3) 


where  m 


■  A  p  u 

e  e  e 


diyat  -  ;B  (4) 

i  -  5  »  <s> 

5.  tot  ■  *  *.  * \/2 


e  -  cy  0 


*TMb  report  uses  the  rotation  of  Reference  j;  definition*  of  eymbole 
are  found  in  LIST  OF  SYMBOLS. 
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Idealized  Geometry  for  Rifle  Gas  System 


1  atm » 


For  computations  thus  far  performed,  8  ■  300°K,  p  ■ 

c  c 

*  0,  v  *  0  at  t  ■  0. 

B  B 

The  quantities  at  the  entrance  to  the  chamber  (e)  remain  to  be  deter¬ 
mined.  This  determination,  however,  entails  the  solution  of  partial 
differential  equations  for  the  unsteady  viscous  heat-conducting  flow 
through  the  duct,  with  initial  and  boundary  conditions  that  require 
matching  this  flow  on  the  one  side  to  conditions  in  the  barrel  (via  the 
port) ,  and  on  the  other  side  to  conditions  in  the  chamber.  The  pro¬ 
cedures  for  obtaining  solution  are  described  in  Reference  1. 

Of  prime  significance  is  the  friction  factor  f  occuring  in  the 
momentum  equation  (also,  by  Reynolds  analogy,  in  the  energy  equation) : 

3u/3t  +  u  3u/3x  +  (1/p)  3p/3x  +  (2  f  u  /D)  *  0 

It  includes  losses  due  to  bends  as  well  as  to  friction  in  the  duct  and 
is  given  by 


f  -  (D/(4l)]  (Ky  +  Ky)  (C) 

2 

where  Ky  and  are  values  taken  from  handbooks.  Since  the  value  of  f 
is  an  estimate  at  best,  it  is  taken  to  be  a  constant,  adjusted  within 
limits  to  give  best  agreement  with  experiment. 

The  analysis  assumes  an  instantaneous  pressure  rise  in  the  barrel 
at  the  port  and  predicts  the  delay  between  this  instant  and  the 
beginning  of  pressure  rise  in  the  chamber. 

111.  STUDY  OF  M-16  GAS  SYSTEM 
A.  Introductory  Remarks. 

Since  all  rifle  gas  systems  can  be  treated  as  the  basic  configur¬ 
ation  of  Figure  1  for  early  times,  it  is  expected  that  sane  general 
Information  could  be  extracted  from  the  study  of  a  particular  system — in 
this  case,  that  of  the  M-16  rifle,  shown  in  Figure  2.  The  bolt  carrier 
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C  (Bolt  Carrier) 
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assumes  the  role  of  the  piston  here#  and  the  annular  bolt  carrier  cavity 
is  the  chamber  of  Figure  1.  Properties  of  the  powder  gas,  dimensions  of 
the  gas  system,  and  loss  coefficients  are  given  in  APPENDIX  A.  Input 
data  for  the  numerical  study  are  drawn  mainly  from  two  rounds  fired  with 
different  ammunition,  which  will  be  referred  to  as  "lot  #1"  and  "lot  #2". 

B.  Input  Data 

The  accuracy  of  the  results  of  the  computation  depends  on  the 
accuracy  with  which  the  input  pressure  is  measured.  In  the  experiments 
static  pressures  were  recorded  on  calibrated  oscillograms  (as  traced  in 
Figure  3).  Unfortunately,  because  of  certain  experimental  factors,  e.g., 
electronic  circuit  delays,  it  has  been  impossible  so  far  to  determine 
accurately  the  initial  pressure  behavior  from  Figure  3.  The  dip  in 
pressure  near  the  peak,  and  possibly  the  peak  dip,  was  caused  by  a 
change  in  ground  potential  of  electrical  equipment.  It  is  necessary  to 
smooth  this  dip  and  other  spurious  signals. 

Computations  shown  in  Figure  4a  indicate  that  changing  the  peak 
pressure  by  13*  has  little  effect  on  the  cavity  pressure. 

Reference  time  (t-0)  in  the  computation  is  taken  as  the  time  at 

which  the  pressure  pulse  appears  at  the  port.  The  computation  assumes 

that  the  pressure  rises  suddenly  to  p  at  the  instant  the  bullet 

niftx 

passes  the  port.  However,  this  time  is  not  well  defined  in  the  experi¬ 
ment.  Between  rise  of  pressure  and  maximum  of  the  pressure  there  is  a 
time  interval  of  about  0.1  x  10  ^  sec.  The  origin  is  thus  uncertain 
by  | 0.05  x  10  | sec.  Experimentally  this  uncertainty  may  be  caused  by 

finite  response  times,  but  also  it  is  noted  that  the  shock  wave  which 
precedes  the  bullet  and  gas  blown  by  the  bullet  cause  an  earlier  rise 
in  pressure. 

To  study  the  effect  of  location  of  reference  time  on  cavity 

pressure,  the  origin  of  time  was  shifted  in  the  computations  by  various 

values  At.  from  the  instant  of  pressure  rise  on  the  port  pressure 
a 

oscillogram.  Results  are  shown  in  Figure  4b.  Best  agreement  occurs  for 
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*  0.1  x  10  sec.  When  this  value  is  applied  as  a  correction  to  the 

measured  delay  between  the  rise  times  of  port  pressure  and  cavity  pres- 
-3  -3 

sure  (*  0.3  x  IQ  sec.),  a  delay  of  0.2  x  10  sec.  results.  As  stated 
in  Section  II,  the  gas  system  model  predicts  the  delay  time.  In  the 

-3 

present  instance  it  predicts  a  delay  of  0.2  x  10  sec.,  in  agreement 
with  the  value  determined  above. 

The  influence  of  difference  in  ammunition  is  indicated  in  Figure  3 
for  two  rounds.  The  greatest  difference  in  pressure  appears  in  the  first 
half-millisecond  of  the  cycle;  the  ultimate  effect  on  cavity  pressure 
will  be  shown  in  Figure  5. 

The  temperature  at  the  port  is  also  required  input  for  the  cal¬ 
culation,  and  is  computed  according  to 

8  (t)/0 (0)  -  [p(t)/£(0)]  (Y  "  1)/Y 

The  initial  temperature  8^(0)  was  not  accurately  determined  for  the 
rounds  fired,  and  an  estimate  was  required.  Fortunately,  however,  cal¬ 
culations  shown  in  Figure  4a  indicate  that  within  limits  (1300*-1500*K  in 
this  .case)  an  exact  knowledge  of  the  temperature  is  not  extremely  crucial. 

C.  Cavity  Pressure 

1.  Comparison  with  Experiment.  Pressure  is  the  variable  of  prime 
interest  here,  and  it  will  be  examined  in  most  detail.  Pressure  his¬ 
tories  computed  for  the  bolt  chamber  are  shown  in  Figure  5.  Agreement 
with  experiment  is  quite  good.  Maximum  pressure  is  reached  in  the  chamber 

-3 

about  1  x  10  sec.  after  the  bullet  has  passed  the  port  station  in  the 
barrel. 

2.  Effect  of  Ammunition.  The  two  input  curves  in  Figure  3, 
representing  the  effect  of  change  of  ammunition,  yield  the  two  theoreti¬ 
cal  cavity  pressure-time  curves  in  Figure  5.  The  difference  in  peak 
values  is  about  fifteen  per  cent. 

3.  Effect  of  Friction  in  Duct.  The  frictional  resistance  to  flow 
in  the  duct  causes  pressure  losses  which  reduce  the  pressure  in  the  bolt 
chamber.  This  resistance  is  incorporated  into  the  model  by  means  of  the 
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"friction  factor"  f_.  Powder  gas  particles  accumulated  from  successive 
firings  adhere  to  the  wall  of  the  duct  and  thereby  produce  progressive 
change  in  the  friction  factor.  The  quantity  f  is  difficult  to  measure 
accurately  and  {as  stated  in  Section  II)  it  must  be  eventually  estimated 
even  after  an  appropriate  value  range  is  determined  with  the  aid  of 
handbook  data. 

In  Figure  6  are  shown  cavity  pressure  histories  for  three  values  of 

f (=.0210,  .0241,  .0270,  where  .0241  is  the  nominal  handbook  value).*  The 
conclusion  is  that  the  pressure  is  not  a  sensitive  function  of  f  in  the 
range  which  yields  optimum  agreement  with  measurements.  On  the  scale  to 
which  the  curves  are  plotted  the  only  notable  differences  can  be  found 
near  the  peak,  where  a  twenty-five  per  cent  variation  in  f  produces  only 

a  seven  per  cent  variation  in  maximum  pressure.  All  these  curves  fall 
below  experimental  pressure  values  on  the  right  hand  side. 

The  f=0  curve  (drawn  to  half  scale  in  Figure  6)  demonstrates  the 
decided  effect  that  the  presence  of  friction  and  heat  conduction  has  on 
the  pressure  in  the  chamber.  Without  inclusion  of  friction  the  pressure 
rises  to  a  maximum  about  two  and  a  half  times  that  of  the  measured  peak 
in  about  two-thirds  of  the  time  taken  to  attain  the  latter,  and  then 
drops  at  a  rate  several  times  that  of  the  measured  pressure.  The  f  ■  0 

curve  can  be  viewed  as  the  limiting  case  for  diminishing  friction  in 
the  duct  flow,  it  is  then  seen  that  the  sensitivity  of  the  pressure  to 
changes  in  f  will  increase  as  f  decreases. 

4.  Frictional  Resistance  to  Piston  Motion.  Inclusion  of  forces 
resistant  to  the  motion  of  the  bolt  carrier,  #  in  Eq.  (3) ,  should 
improve  agreement  between  theory  and  experiment .  An  estimate  of  the 
magnitude  of  the  forces  involved  can  be  gained  by  a  knowledge  of  the 
discrepancies,  as  in  Figure  6.  However,  the  friction  force  law  is  not 
presently  known  in  this  situation. 

A  simplified  friction  law  was  assumed  and  two  cases  computed  with 

iti  namely,  that  the  piston  does  not  move  until  the  cavity  pressure 

reaches  a  preassigned  value  p  ,  then  it  moves  against  a  force 

s 


^Because  the  friotion  factor  depends  on  the  particular  gae  ay  atm, 
through  Bq .  (8),  the  value  of  f  can  vary  widely  from  one  gun  eyetm  to 
another. 
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Figure  6.  effect  of  Friction  Factor  f  on 

Cavity  Fraaaura  (M-16,  Lot  12  Mound) 


4  =  -  Ap^A^,  where  0  <  A  ■  const.  <  1.  Results  shown  in  Figure  7 
indicate  that  this  model  does  not  uniformly  improve  agreement  with 
experiment;  there  is  a  tendency  to  produce  peak  pressures  too  high  in 
the  process  of  closing  the  gaps  on  the  sides  of  the  curves.  More 
study  will  be  required  on  this  subject, 

5.  Effect  of  Jamming.  To  determine  the  effect  of  gas  system 

jamming  on  cavity  pressure,  a  computation  was  performed  in  which  v  was 
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set  equal  to  zero  for  t  1  2.05  x  10  sec.  The  result  is  shown  in 
Figure  8  for  a  particular  round.  A  pressure  rise  occurs  immediately 
after  the  piston  stops,  then  the  pressure  decreases  at  a  slower  rate 
than  for  the  unimpeded  piston.  Qualitatively  this  behavior  has  been 
confirmed  by  experiments. 

D.  Cavity  Temperature  and  Density 

Temperature  and  density  histories  in  the  cavity  are  plotted  in 
Figures  9  and  10.  It  is  seen  that  the  temperature  reaches  peak  value 
very  rapidly;  the  cooling  which  then  follows  is  due  to  the  expansion 
of  the  gases  in  the  expanding  cavity  as  well  as  to  the  decrease  of  IS 

9 

with  time.  This  cooling  brings  the  temperature  down  to  the  range  of 
its  initial  value  (in  this  case,  below  it)  within  the  first  2.4  x  10~3 
sec. 

The  density  attains  its  maximum  at  approximately  1  x  10~3  sec., 
and  at  2.4  x  10  3  sec.  it  is  still  about  twenty  times  its  initial  value. 

The  rates  of  fall  of  temperature  and  density  are  increased 
slightly  when  the  vent  is  uncovered. 

A  survey  of  the  calculations  shows  that  temperature  and  density 
are  generally  affected  less  than  pressure  by  the  variations  of  the 
parameters  discussed  above  in  Subsection  C.  The  limiting  case  of  zero 
heat  conduction  and  friction  (f  ■  0)  again  has  a  marked  effect;  the 
peak  temperature  Increases  by  a  factor  of  four,  and  the  peak  density 
is  cut  by  about  one- half  of  its  value. 
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Figure  8.  Effect  of 
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Figuro  9.  TMgmtuw  la  Cavity  (H-16,  Lot  12  Jtoond) 
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E.  Motion  of  Piston 


The  motion  of  the  bolt  carrier  relative  to  the  bolt  is  depicted 
in  Figure  11,  where  relative  velocity  and  bolt  carrier  displacement 

_3 

are  shown.  In  about  2.4  x  10  sec.  the  piston  has  moved  .0076  m.  and 
attained  a  relative  velocity  of  6.3  m/sec. 

The  velocity,  v  ,  of  the  bolt  carrier  at  time  t_,  when  it  collides 
BE 

with  the  bolt,  is  of  interest  here,  since  it  is  a  measure  of  the  momentum 
applied  to  the  reloading  system.  Variation  of  the  friction  factor  f 
from  0.021  to  0.027  produces  less  than  three  percent  change  from  the  t 

E 

and  v  of  Figure  11.  The  v  for  the  Lot  #1  ammunition  round  differs  by 

B£  DL 

less  than  six  percent  from  that  for  the  Lot  #2  round. 

Frictional  resistance  produces  retardation  ;  for  the  case  p  =  0, 
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A  =  0.3  it  takes  2.8  x  10  sec.  for  the  piston  to  attain  a  velocity  of 
5.0  m/sec. 

With  friction  and  heat  conduction  absent  from  the  duct  flow  (f  =  0) 

_3 

a  piston  velocity  of  8.73  m/sec.  is  reached  in  1.7  x  10  sec. 

IV.  AR-18  GAS  SYSTEM 

A.  Discussion 

The  AR-18  rifle  is  basically  the  same  gun  as  the  M-16  rifle.  It 
has  the  same  barrel  and  location  of  port  as  the  latter?  however,  the 
gas  system  is  constructed  differently.  A  sketch  of  the  AR-18  gas 
system  mechanism  is  presented  in  Figure  12  at  two  instants  in  its 
cycle:  i)  beginning  of  the  piston  motion,  and  ii)  uncovering  of  the 
vents.  The  dimensions  of  the  system  and  the  loss  coefficients  are 
given  in  APPENDIX  A. 

A  numerical  study  of  the  AR-18  gas  system  was  carried  out,  using 
input  data  employed  in  the  M-16  study  because  of  the  basic  similarity 
of  the  two  rifles. 


30 


Velocity  and  Displacement  History  of  Piston  (M-16,  Lot  #2  Round  ) 
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Figure  12.  Sketch  of  AR-18  Rifle  Gee  Syeten 
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B.  Calculations 


The  computational  results  for  a  round  fired  in  the  AR-18  rifle  are 
shown  in  Figures  13  and  14.  Qualitatively  the  variables  behave  in  the 
same  manner  as  the  corresponding  ones  in  the  M-16  rifle.  A  comparison 
with  measurements  of  displacement  vs  time  and  final  velocity  of  the 
piston  (the  only  experimental  data  available)  shows  good  agreement 
between  theory  and  experiment. 

C.  Comparison  with  M-16 

While  the  dynamic  and  thermodynamic  variables  are  all  of  the  same 
order  of  magnitude  for  the  M-16  and  AR-18  gas  systems,  there  are 
quantitative  differences  due  to  differences  in  size  of  corresponding 
parts . 

Because  of  the  shorter  duct  length,  the  delay  time  (between 
pressure  rise  in  pert  and  that  in  cavity)  of  the  AR-18  is  about  one- 
fifth  of  that  of  the  M-16.  Venting  produces  a  steeper  drop  in  the 
thermodynamic  quantities  in  the  AR-18,  which  has  a  venting  area  about 
four  times  that  of  the  M-16.  The  peak  pressure  of  the  AR-18  is  about 
one  and  a  half  times  that  of  the  M-16,  and  the  peak  temperature  is 
about  100  °K  higher  in  the  AR-18. 
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APPENDIX  A:  TABLES  OF  PARAMETERS* 


Table  A-I.  Properties  of  Powder  Gas** 

=  1.74  x  103  (m/sec) 2/°K 

cy  -  1.38  x  103  (m/sec) 2/°K 

ic  -  0.831  x  10-1  (kg  m)/(sec3  °K) 

R  *  0.400  x  103  (m/sec) 2/°K 
gas 

Y  -  1.26 

?  -  1.24 

-5 

y  *  4.80  x  10  kg/ (m  sec) 


Table  A-2.  Dimensions  of  M-16  Gas  System 


-  1.269  x  10"4  »2 

-  0.6605  x  10~5  m2 

-  0.4382  x  10"  5  m2 

-  0.6605  x  10‘5  m2 

-  0.6533  x  10'5  *2 

-  0.2900  x  10’2  m 
*  0.3600  ■ 


r 

v 


*BE 

*Bv 

(C/D) 


-  0.4366  kg 

•  0.1442  x  10"2  m 

-  0.7600  x  I0"6  m3 

-  0.7600  x  10"2  m 

-  0.6158  x  10’2  m 
■  0.015 


'Values  of  ponder  gat  properties  and  dimtntiont  of  gat  ay  stem  supplied 
by  Interior  Ballietiee  Laboratory. 

"Estimates  bated  on  the  interior  ballietioe  oosputation  of  the 
composition  of  the  ponder  gat. 
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Table  A-III.  Loss  and  Flow  Coefficients  of  M-16  Gas  System 


a  =  0.62 
e 

Kb  *  5.04 

a  =  0.64 

P 

Kp  =  6.95 

a  =  0.62 

g 

f  =  0.0241 

Table  A- IV.  Dimensions 

of  AR-18  Gas 

System 

A 

c 

- 

0.732  x  10~4  sn2 

0.463  kg 

A 

e 

= 

0.766  x  10~5  m2 

r  = 

V 

0.299  x  10-2  m 

A  . 
min 

= 

0.195  x  10~5  m2 

V  . 

Cl 

0.690  x  lO-6  m' 

A 

P 

= 

0.766  x  10~5  m2 

XBE 

1.57  x  10“2  m 

A 

V 

= 

2.80  x  10“5  m2 

XBv 

0.970  x  10“2  m 

5 

= 

0.310  x  10“2  m 

(e/D)  “ 

0.015 

I 

— 

0.0660  m 

Table  A-V.  Loss  and  Flow  Coefficients  of  AR-18  Gas  System 

ag  *  0.62  Kg  -  3.248 

ap  *  0.64  Kg  *  1.1815 

a  -  0.62 

g 


f  -  0.0525 


APPENDIX  B:  CONVERSION  TABLE 


Area: 

1 

m2  =  10.764  ft2  =  1550  in2 

Density: 

1 

kg/m3  =  0.06243  lb/ft3 

Energy: 

1 

kg  (m/sec)2  =  0.94787  x  10~3  BTU 

Force : 

1 

N  =  1  kg  m/sec2  =  0.22482  lbforce 

Length : 

1 

m  =  3.281  ft  =  39.37  in. 

Mass: 

1 

kg  =  2.205  lb 

mass 

Pressure: 

1 

N/m2  =  1  kg/ (m  sec2)  -  0.14511  x  10' 

Velocity: 

1 

m/sec  =  3.281  ft/sec  *  39.37  in/sec 

Volume : 

1 

m3  =  35.315  ft3  *  61.024  x  103  in3 

6  2 

1  atm.  press.  «*  0.1013  x  10  kg/(m  secz)  *  14.69  psi 
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